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INTRODUCTION 


This  report  contains  a  description  of  an  experimental  videotape  reading  sys¬ 
tem  developed  at  the  White  Sands  Missile  Range  Instrumentation  Directorate 
computer  lab  for  the  investigation  of  image  processing  and  pattern  recogni¬ 
tion  concepts.  The  VRS  is  presently  being  used  to  study  the  concept  of  accur¬ 
ately  determining  the  aspect  angles  of  a  target  from  one  frame  of  video.  The 
ability  of  accurately  finding  the  position  of  a  target  from  one  frame  of  video 
Is  useful  In  extracting  a  data  product  from  a  videotape  when  there  is  tape 
available  from  only  one  station.  Such  a  system,  made  into  a  real-time  hard¬ 
ware  machine,  would  also  have  applications  in  fire  control  of  high-energy  lasers, 
since  the  aiming  of  such  devices  requires  that  exact  knowledge  about  the  posi¬ 
tion  of  the  target  be  available  so  that  energy  can  be  deposited  at  a  critical 
point  of  the  target.  This  experimental  system  is  useful  as  a  test  bed  for  con¬ 
cepts  that  will  have  applications  both  in  extracting  a  data  product  to  be  used 
by  customers  of  WSMR  and  as  a  model  for  a  hardware  machine  that  would  be  used 
both  for  real-time  tracking  and  for  fire  control  of  new  weapons  technology. 

SYSTEM  DESCRIPTION 

The  data  flow  of  this  system  is  as  follows: 

‘  A  video  tape  of  a  mission  is  taken  at  a  station. 

*  The  frames  to  be  read  are  put  on  a  video  disk  which  is  attached  to 
an  image  analyzing  system  capable  of  digitizing  the  frames  in  the  video  disc. 

’  These  frames  are  digitized  and  put  into  data  files  with  another 
file  containing  all  the  file  names  of  the  frames  which  are  of  interest. 

The  software  then  processes  the  data  in  the  following  sequence: 

1.  Read  in  file  containing  names  of  files  to  be  processed. 

2.  Read  in  first  file  and  do  preprocessing  on  it  until  completely  done. 


begin 

repeat 

cobegin 

begin 


3a.  Read  next  file  and  do  preprocessing  on  it. 


3b.  Make  a  contour  of  the  previously  processed  file  and  do  the  classi¬ 
fication. 

end; 

end; 

until  eof; 


4.  Finish  off  classifying  last  one  read  in 


Processing  begins  by  first  doing  pixel  level  operations.  The  classifica¬ 
tion  is  done  by  making,  a  line  drawing  of  the  plane  or  rocket  to  be  analyzed, 
and  comparing  it  against  a  previously  stored  line  drawing  library  made  from 
views  of  the  object  in  question  at  different  angles.  Before  a  contour  of 
the  target  (Fig.  la)  can  be  made,  points  which  are  possible  candidates  for 
edge  points  must  be  identified.  Since,  typically,  scenes  that  we  process 
are  very  noisy,  we  begin  by  doing  a  three-by- three  averaging  to  every  point 
in  the  scene.  After  this,  a  moment  edge  detector  is  used  to  assign  to  each 
point  in  a  scene  a  value  which  reflects  the  probability  that  a  point  is  an 
edge  point  (Fig.  lb).  A  threshold  is  chosen  by  the  operator  and  all  points 
classified  as  possible  target  points  are  assigned  a  zero  and  all  others  a 
one.  The  computer  then  makes  and  displays  a  contour  of  the  entire  scene 
with  different  polygonal  segments  being  assigned  different  values  (Fig.  2). 
The  operator  chooses  the  number  of  segments  which  make  up  the  target,  and 
the  computer  writes  the  segments  out  in  a  file.  This  file  is  then  modified 
by  the  use  of  interactive  graphics  programs  (Figs.  4  and  5).  The  result 
(Fig.  lc)  is  compared  against  the  library  of  stored  views  (Fig.  3),  the 
best  match  is  found  (Fig.  6)  and  the  angular  data  needed  is  read  from  the 
coordinate  system.  A  description  of  the  operations  that  take  place  is  thus: 

begin 

1.  Read  in  file  containing  names  of  files  to  be  processed. 

2.  Read  in  first  file  and  do  preprocessing  on  it  until  completely  done. 

repeat 

cobegin 

begin 

3A.  Read  next  file  and  do  preprocessing  on  it. 

end; 

begin 

3B.  Make  a  contour  by  the  following  process: 

a.  Using  a  histogram,  computer  chooses  a  threshold  for  the 
moment  file  of  the  original,  and  displays  a  contour  based  on  this  threshold. 

b.  Is  this  contour  acceptable? 

c.  while  contour  not  acceptable  do 
begin 

♦Obtain  new  threshold  from  operator. 

♦Draw  contour 

♦Is  contour  acceptable? 

end 

d.  Let  operator  choose  segments  that  will  be  used  to  construct 

target. 

e.  Display  segments  chosen  by  the  operator  and  modify  them  as  the 
operator  instructs. 

f.  Calculate  the  Fourier  descriptors,  normalize  and  do  classifi¬ 
cation. 

enu? 

until  eof; 

4,  Finish  off  classifying  last  one  read  in 


J 


end. 
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Before  a  contour  of  the  scene  can  be  made,  points  which  are  possible  candi¬ 
dates  for  edge  points  must  be  identified.  Since,  typically,  the  scenes  that 
we  process  are  very  noisy  we  begin  by  doing  a  three-by-three,  averaging  to 
every  point  in  the  scene.  After  this  a  moment  edge  detector  is  used  to 
assign  to  each  point  in  a  scene,  a  value  which  reflects  the  probability  that 
a  point  is  an  edge  point.  As  shown  in  Appendix  A  this  sequence  of  steps 
increases  the  probability  of  detecting  edge  points.  The  next  step  is  to 
do  another  averaging  operation  on  the  moment  file  with  the  purpose  of  in¬ 
creasing  the  connectedness  of  the  edge  points. 


The  software  to  accomplish  the  preprocessing  was  written  with  two  ends  in 
mind;  one  was  that  this  software  would  be  a  model  for  a  hardware  module  to 
be  built  later,  and  the  other  was  that  execution  time  be  reduced  by  overlap¬ 
ping  input/output  with  processing.  Figure  7  illustrates  how  the  software  is 
set  up.  The  programs  read  l,  write  i  and  process  are  passive  programs  in 
that  they  suspend  themselves  immediately  after  doing  some  initialization  oper¬ 
ations.  These  consist  of  bookkeeping  operations  such  as  setting  input  file 
name,  output  file  name,  and  setting  up  parameters  so  that  the  proper  buffer 
is  accessed  each  time  a  program  is  activated.  The  program  which  drives  these 
passive  programs  is  called  main  l.  It  runs  the  needed  programs  and  synchro¬ 
nizes  them  via  the  use  of  global  event  flags.  After  the  preprocessing  is 
finished  it  initiates  the  next  step  in  processing  by  its  call  to  arrows. 


A  typical  frame  is  processed  by  main  I  in  the  following  way:  First  the  pro¬ 
grams  READi ,  writei  and  process  are  loaded  into  memory.  They  do  whatever 
initialization  is  necessary  and  then  suspend  themselves.  There  are  two  in¬ 
put  buffers  that  will  be  used  by  readI  to  store  the  data  to  be  processed, 
and  two  output  buffers  where  processed  data  is  put  and  from  where  the  program 
writei  writes  the  data  out  onto  the  disk,  maini  first  has  the  two  input  buf¬ 
fers  (lines  15  -  IB)  filled  by  the  two  activations  of  readi  done  by  two  calls 
to  resume (readi) .  readi  automatically  processes  the  buffers  in  an  alternate 
manner  as  do  proces  and  writei.  The  buffers  are  initially  set  up  (lines  21  - 
24)  so  that  the  remalnina  processing  can  be  done  concurrently  (lines  25  -  35). 
In  the  do  loop  there  are  waits  for  flags  to  be  set  that  indicate  that  each  of 
the  programs  Involved  are  finished.  The  rest  of  maini  finishes  up  with  the 
buffers  that  need  to  be  processed  and  written  out.  On  line  45  it  starts  the 
next  step  for  this  frame  by  its  call  to  arrows. 


0015 

0016 

0017 

0018 

0019 

0020 

C 

0021 

0022 

0023 

C 

c 

c 

0024 


CALL  RESUME  (READI)  FILL  IN  BUFFER  #2 

CALL  WAITFRC36) 

CALL  CLREF (36) 

CALL  RESUME (READI) 

CALL  WAITFR  (3f>) 

CALL  CLREF  (36)  (  BUFFERS  #1  and  #2  FULL 

CALL  RESUME (PR0CE1 , ) 

CALL  WAITFR ( 37) 

CALL  CLREF (37) 

AT  THIS  POINT  INBUF#1  AND  INBUF#2  ARE  FILLED 
AND  #2  HAS  BEEN  COPIED  OVER  TO  OUTBUF  #2 

CALL  WAITFR (42) 
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0025 

DO  30  L-1,2 

0026 

CALL 

RESUME (READI) 

0027 

CALL 

RESUME (WRITE1) 

0028 

CALL 

RESUME  0  ROCEl ) 

0029 

CALL 

WAITFR ( 38 ) 

003'' 

CALL 

WAITFR { 37) 

0031 

CALL 

WAITFR (36) 

0032 

CALL 

CLREF (38) 

0033 

CALL 

CLREF (37) 

0  34 

CALL 

CLREF (36) 

0035 

30 

CONTINUE 

0036 

CALL 

RESUME { WRITE 1) 

00  3  7 

CALL 

RESUME (PROCE1) 

0038 

CALL 

WAITFR (3o) 

0039 

CALL 

W/  ITFRC37) 

0040 

CALL 

CLREF (3M) 

0041 

CALL 

C!  REF  ;.J7) 

0042 

CALL 

SETEF (40) 

0043 

CALL 

RESUME (WRITEl) 

0044 

CALL 

WAITFR ( 38) 

0045 

CALL 

REQUES (ARROWS) 

read i  has  access  to  a  conrnon  global  area  (line  8)  where  It  reads  In  the 
data  to  be  processed.  After  initializing  (lines  19  -  20)  it  suspends  It¬ 
self  until  activated.  When  activated  it  places  the  last  six  rows  it  has 
read  in  the  first  six  positions  of  the  buffer  (lines  25  -  27)  to  be  pro¬ 
cessed.  It  then  reads  32  rows  and  stores  the  last  six  into  store.  The 
bookeepinn  for  the  change  of  buffers  is  done  (lines  32  -  34).  maini  is 
signalled  that  readi  Is  finished  by  the  call  to  setef  and  there  is  a  jump 
to  21  which  suspends  readi. 


JoOl 

I  RC  .PAM  READI 

0002 

INTEGER* 2  SWITCH, POINT, WAKEUP  .FINISH  , OFFSET 

0003 

INTEGER*2  STORE (-383:0) 

0004 

INTEGER* 2  INTIN (-383:2048, 2) , INTOUT ( -383 : 2048 , 2 

000  5 

LOG I CAL *1  INEUFF (128,-5:32,2)  l  INPUT  BUFFER 

0006 

LOGICAL* 1  OUBUFF (128,-5:32,2)  1  OUTPUT  BUFFER 

0007 

LOGICAL  *1  BlQ (25) 

0003 

COM'IO!.  /DTA/INBUFF  .OUBUFF 

0009 

EQUIVALENCE  ( INTIN (- 383 , 1 ) , INBUFF ( 1 ,- 5 , 1 ) > 

0010 

EQUIVALENCE  (OUBUFF ( 1 , -5 , 1 ) , INTOUT ( -3S 3 , 1 )) 

0011 

DATA  WAKEUF , FINISH/3 3, 36/ 

0019 

POINT-1 

0020 

SWITCH- 2 

0021 

CALL  CLREF (40) 

0022 

CALL  SETEF (41) 

0023 

20 

CALL  SUSPND 

0024 

READ (lO1 POINT) ( INTIN ( J , SWITCH) ,J-1,2048) 

0025 

DO  15  I--383 ,0 

0026 

INTIN (I .SWITCH)-  STORE ( I ) 

0027 

15 

CONTINUE 
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0028 

DO  16  1—383,0 

0029 

STORE ( I ) -  INTIN (  2048+1, SWITCH) 

0030 

16 

CONTINUE 

0031 

30 

CONTINUE 

0032 

INTER-SWITCH 

0033 

IF ( INTER. EQ,1)  SWITCH- 2 

0034 

IF(INTER.EQ.2)  SWITCH-1 

0035 

IF (POINT, EQ. 5)  GOTO  40 

0036 

CALL  SETEF (FINISH) 

0037 

GOTO  20  I  GO  WAIT  TILL  AWOKEN 

0038 

40 

CALL  SETEF (39) 

0039 

CALL  SETEF (FINISH) 

0040 

END 

The  next  step  after  the  files  have  been  read  Is  to  do  the  averaging  and 
edge  detection.  Again  proces  has  access  to  the  global  common  area  DTA. 

It  initializes  Itself  and  then  suspends  Itself  and  waits  for  maini  to 
activate  It  when  needed.  The  processing  (lines  13,  14,  15)  consists  of 
an  averaging  operation,  an  edge  detection  (moment)  and  another  averaging. 
The  processing  Is  done  from  inbuf  to  outbuf  (AV6),  outbuf  to  inbuf  (mo¬ 
ment),  and  then  inbuf  to  outbuf.  The  bookkeeping  to  allow  alternate  buf¬ 
fers  to  be  processed  Is  then  done;  the  program  suspends  Itself  and  then 
waits  for  the  next  call. 

OOOl  PROGRAM  PROCESS 

0002  INTEGER*2  SWITCH,POINT,WAKEUP,AVERAG,FINAVG 

0003  REALM  Ml, M2, MX, MV 

0004  L0GICAL*1  INBUF(l28,-5:32,2)  l  INPUT  BUFFERS 

0005  L0GICAL*1  OUTBUF  (128,-5:32,2)  .'  OUTPUT  BUFFERS 

0006  COMMON  /DTA/  INBUF, OUTBUF 

0007  DATA  WAXEUP .FINAVG/34, 37/ 

0008  SWITCH  =  2 

0009  KDUNT=0 

0010  CALL  CL REF (WAKEUP) 

0011  20  CALL  SUSPND 

0012  ROUNT-KOUNT+1 

0013  CALL  AVERAG (SWITCH, 31) 

0014  CALL  MOMENT (SWITCH, 29) 

0015  CALL  AVERAG (SWITCH, 27) 

0016  CALL  SETEF (FINAUG) 

0017  INTER- SWITCH 

0018  IF ( INTER .  EQ .  1 )  SWITCH- 2 

0019  IF(INTER.EQ.2)  SWITCH-1 

0020  IF (KOUNT . EQ , 4)  GOTO  21 

0021  GOTO  20 

0022  21  CONTINUE 

0023  END 


0001 

SURROUTINE  AVERAG (SWITCH, L) 

0002 

L0GICAL*1  INBUF (128,-5:32,2) ,OUTBUF (128 ,-5 : 32 , 2) 

0003 

INTEGER* 2  Rl,R2,R3 

0004 

INTEGER*2  SWITCH, L 

0005 

COMMON  /DTA/  INBUF, OUTBUF 

0006 

DO  10  J=-4 , L 

0007 

D'  20  1=1,127 

0008 

R1=R2 

0009 

R2-R3 

0010 

10  =  INBUF (1+1 , J-l , SWITCH) .AND  255 

0011 

11  =  INBUF (I+1,J, SWITCH) .AND. 255 

0012 

12  =  INBUF (I+l.J+l  SWITCH) .AND. 255 

0013 

R3  =  (IQ+I1+I2) /3 

0014 

IAVG=  (RI+R2+R3J/3 

0015 

OUTBUF  1,1,  J-l,  SWITCH)  =  IAVG .  AND .  2  55 

0016 

20 

CONTINUE 

0017 

10 

CONTINUE 

0018 

RETURN 

0019 

END 

0001 

SUBROUTINE  MOMENT (SWITCH , L) 

0002 

L0GICAL*1  INBUF (128, -5: 32, 2) , OUTBUF (128 , -5 : 32 , 2) 

0003 

INTEGER* 2  SWITCH 

0004 

COMMON  /DTA/  INBUF, OUTBUF 

0005 

DO  30  J=-4,L 

0006 

DO  40  1=2,127 

0007 

10- 1 3 

0008 

11=14 

0009 

12=15 

0010 

13=16 

0011 

14=17 

00  12 

15=  lb 

0013 

I6=OUTBUF(I+l, J-l, SWITCH) .AND. 255 

0014 

I 7-OUTBUF (I+1,J, SWITCH) .AND. 255 

0015 

I8=OUTBUF (I+1,J+1, SWITCH) .AND. 255 

0016 

XM=FLOAT (5* (10-18) +4* (I1+I3-I5-I7) ) 

0017 

YM=FLOAT (5* (I6-I2)+4* (I3+I7-I1-I5) ) 

0018 

M=SQRT (XM**2+YM**2) 

0019 

INBUF (I, J-l , SWITCH) =M. AND.  255 

0020 

40 

CONTINUE 

0021 

30 

CONTINUE 

0022 

RETURN 

0023 

END 

The  averaqe  that  Is  done  Is  an  unweighted  average.  The  edge  detector  used 
is  a  moment  operator  which  has  been  shown  to  perform  well  In  the  presence 
of  noise.  The  next  program  that  Is  called  Is  writei.  The  data  structure 
here  are  the  same  as  those  used  for  readi  with  the  same  global  common  area 
being  used.  It  also  suspends  Itself  and  waits  to  be  activated. 


0001 

PROGRAM  WRITEl 

0002 

INTEGER* 2  SWITCH, POINT, WKEUP,FINWRI , OFFSET 

0003 

INTEGER* 2  INTOUT  (-383 -.2048, 2) 

0004 

LOGICAL* 1  INBUF (128,-5:32,2)  ! 

INPUT  BUFFERS 

0005 

L0GICAL*1  OUTBUF (128  -5:32,2)  . 

OUTPUT  BUFFERS 

0006 

L0GICAL*1  BIO (25),  CHARAC 

0007 

COMMON  /DTA/ INBUF, OUTBUF 

0008 

EQUIVALENCE  (INTOUT (-383 , 1 ) , OUTBUF (1,-5, 1) ) 

0009 

DATA  WAKEUP,  FINWRI/35 , 38/ 

0018 

POINT-1 

0019 

c 

SWITCH-2 

0020 

V 

CALL  CLREF (WAKEUP) 

0021 

CALL  SETEF (42) 

0022 

20 

CALL  SUSPND 

0023 

WRITE ( 11' POINT) { INTOUT (J, SWITCH, 

,J=-3B3,1664) 

0024 

30 

p 

CONTINUE 

0025 

W 

INTER- SWITCH 

0026 

IF ( INTER. EQ.l)  SWITCH- 2 

0027 

IF  (INTER. EQ. 2)  SWITCH-1 

0028 

CALL  READEF (40 , LCQDE) 

0029 

IF(LC0DE.EQ.2)  GOTO  46 

0030 

CALL  SETEF (FINWRI) 

0031 

GOTO  20 

0032 

46 

CONTINUE 

0033 

INBUF (1,-5, 1)  »  CHARAC 

0034 

DO  47  1=1,25 

0035 

INBUF (I+l,-5,l)=  B10(I) 

0036 

47 

CONTINUE 

0037 

CALL  SETEF (FINWRI) 

0038 

END 

A  model  for  a  hardware  realization  of  this  software  Is  given  In  Figure  8. 
Mere  each  of  the  circles  would  be  a  cpu  together  with  some  local  memory. 
They  would  be  passive  and  controlled  by  a  cpu, main.  The  squares  would 
correspond  to  buffers  accessed  by  cpu's  as  Indicated.  There  are  standard 
hardware  methods,  such  as  Interrupts  and  flags,  that  can  be  used  for  the 
synchronization  which  Is  done  In  the  software  model. 
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CONSTRUCTION  OF  A  CONTOUR  FROM  THE  MOMENT  FILE 

After  the  moment  file  pic  Is  created  a  threshold  "T"  must  first  be  found 
such  that  a  contour  of  the  target  will  be  included  in  the  set. 


{  PlC(i,j):  i=l..n,J=l..N,PIC<i,j)  <  T  } 


To  find  t  we  first  compute  the  histogram  of  the  moment  file  created.  A 
number  p  needs  to  have  been  chosen  beforehand  which  represents  the  percentage 
of  the  scene  points  which  are  target  points.  The  first  point  for  which 


/  histo  >  P 

A 


is  found  and  used  as  the  value  of  t.  It  has  been  found  that  p=  15.  works 
well  In  cases  where  the  target  is  a  small  part  of  the  scene  and  p=  25.  does 
well  when  the  target  is  a  large  percentage  of  the  scene.  In  Figure  9  there 
are  four  originals  that  will  be  reduced  to  a  contour.  The  result  of  pre¬ 
processing  this  data  is  in  Figure  10.  The  problem  now  is  to  find  a  t  such 
that  the  target  will  be  separated  from  the  background.  If  we  look  at  the 
raw  histograms  (Figure  111)  we  can,  in  some  cases,  guess  at  where  the  thresh¬ 
old  should  be  chosen,  assuming  that  there  is  one  distribution  for  the  target 
and  another  for  the  background.  The  background  distribution  is  centered 
about  the  maximum  of  the  histogram  while  the  target  distribution  is  part  of 
the  tall  of  the  histogram.  Thus  it  is  reasonable  to  suppose  that  the  target 
points  constitute  a  certain  percentage  of  the  points  to  the  right  of  some 
value.  Experiments  have  shown  that  the  proper  value  for  this  percentage  is 
between  15  and  25,  depending  on  the  size  of  the  target.  Figure  11b  is  a 
figure  found  from  11a  by  graphing 


255 

Sum  (X)  =  /  histo. 

1  A  1 


for  each  histogram  of  11(1).  From  this  graph  we  see  that,  as  the  contrast 
decreases,  the  threshold  to  be  chosen  decreases,  a  procedure  that  agrees  with 
our  Intuition.  We  can  also  see  that,  when  the  target  size  is  large,  the 
graph  is  radically  different  than  when  the  target  is  small.  Using  11(11) 
and  p  =  23,  we  obtain  the  contours  of  Figure  12.  The  computer  is  set  to 
threshold  at  p  =  23,  the  contour  appears  on  the  screen;  and  the  operator  can 
reject  this  contour  and  request  a  new  one  based  on  an  operator  supplied  value 
for  p.  One  choice  of  p  does  not  always  produce  closed  contours  of  the  target; 
and  this  Is  why  operator  Intervention  is  required  at  this  point.  As  this 
system  stands  now,  p  Is  set  by  the  operator  on  the  initial  frame  and  used  for 
subsequent  frames  until  the  operator  Intervenes. 
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Figure  10a.  Moment  file  of  FI 02 
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Figure  10b. 


Moment  file  of  FI 02  flyin 
in  front  of  a  mountain. 
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F-.yure  10c.  Moment  file  of  Hawk 
missile. 


Figure  lOd.  Moment  file  of  FI 02 ,  sun¬ 
spots  on  wings  and  nose 


Figure  11.  Raw  histogram  Data  (I)  from  moment  file.  Integrated  Data  (II)  from 
raw  histograms. 
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Figure  12a.  Contour  found  from  Mo¬ 
ment  file  (Fig.  10a.) 

V  * 


Figure  12b.  Contour  found  from  Moment 
file  (Fig.  10b). 


Figure  12c.  Contour  found  from  Mo¬ 
ment  file  (Fig.  10c). 


Figure  12d.  Contour  found  from  Moment 
file  (Fig.  lOd). 
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Once  the  value  of  t  has  been  determined  so  that  the  target  lies  in  the  set 

{PIC(i,  j)  :i=l. .  .n,  j=l.  .m  ,  PIC(i,j)<T  } 

We  are  in  a  position  to  begin  making  the  contours  from  which  the  targets 
will  be  extracted.  Again  since  these  programs  are  intended  to  be  used  as 
models  for  a  future  hardware  implementation  the  processing  follows  a  sequen¬ 
tial  order.  That  is,  the  rows  are  processed  from  l  to  n  with  the  original 
gray  values  used  just  once  in  the  processing.  It  is  commonly  known  that  to 
construct  a  fast  contour  plotting  algorithm  the  Freeman  code  must  be  dropped 
as  soon  as  possible.  We  use  the  Freeman  code  to  record  the  direction  as  we 
go  along  but  convert  to  a  polygonal  representation  as  soon  as  it  is  possible, 
as  when  there  occur  two  consecutive  Freeman  codes  with  the  same  direction. 

This  is  the  first  step  in  the  contour- forming  process  and  results  in  a  set  of 
line  segments  which  are  specified  by  their  endpoints.  The  next  step  is  join¬ 
ing  the  elements  into  contours. 

The  program  that  performs  the  operations  described  above  is  the  program  arrow. 
First  of  all,  six  rows  to  be  processed  are  read  into  ammoun (128, 6) .  One  pair 
of  rows  is  processed  at  a  time  with  the  results  of  the  processing  being  put 
into  stor(128, 2)  (lines  52  -  63).  The  assignment  of  directions  begins  by  — 
first,  thresholding  two  rows  of  ammoun,  with  a  point  being  assigned  a  zero  if 
the  average  of  a  two-by-two  neighborhood  is  greater  than  t  and  a  one  if  it  is 
less  than  t  (lines  59  -  61) .  Beforehand  stor(128,d  is  first  overwritten  by 
the  last  row  processed  which  had  been  put  in  stor(128,2)  (lines  65  -  66). 


OOOl 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 


PROGRAM  ARROW 

INTEGER*2  MAG(384) ,ANG(384) , CENTER, THRSHM,X 
INTEGER* 2  U,V,POINT3,AO,Al,A2 , A3 
INTEGER*2  P0INT1 ,P0INT2 , VAR, SUM, SIGN 
INTEGER*2  hO ,H1 , H2 ,H3 , STOR (128 , 2) ,AVERAG, THRESH 
L0GICAL*1  Y,ANS, INV,FLAG 

LOGICAL*!  AMMOUN (128,6) , ANGLE (128,6) , NEW (128, 5) 
L0GICAL*1  B10(25) ,NAMETE(26) 

REAL *4  XI , SEGMN 

INTEGER* 2  LINES (500, 3 ) COORDI (1500,6) 

INTEGER* 2  POINSI 

INTEGER*2  POINTL,POINTH ,POINTT,COL 
COMMON  INV, FLAG, SIGN 

COMMON  LINES, P0INTL,P0 INTO, PO INTO, COL, 128  ' 

COMMON  /DTA/COORDI 
EQUIVALENCE  (MAG, AMMOUN) 

EQUIVALENCE  (NAMETE, COORDI) 

EQUIVALENCE  (ANG, ANGLE) 

DATA  Y/89/ 

DATA  SEGMN/6RSEGMEN/ 

DATA  IJ/O/ 

DATA  AK/1. 4111764/ 

CALL  ERRSE T (37,. TRUE . , . FALSE .,. FALSE ., . FALSE .,31) 
I=NAMETE (1) 

DO  999  M=l, 25 

BIO (M)  =  NAMETE (M+l) 
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0027 

999 

CONTINUE 

0028 

CALL  ASSIGN (7, BIO, I) 

0r'29 

DO  947  L=1 , 20 

0030 

COORDI (L , 1)  =0 

0031 

COOADI(L,2)=0 

0032 

COORDI (L, 3 )=0 

0033 

COORDI (L, 4 )=0 

0034 

COORDI (L, 5) =0 

0035 

COORDI (L, 6 )=0 

0036 

947 

CONTINUE 

0037 

POINSI= 1 

00  38 

DEFINE  FILE  7 (0, 64 , V, POINT2) 

0039 

POINT2=l 

0040 

CALL  DEVIAT (THRESH, P0INT2) 

0041 

POINT2=l 

0042 

IJ=0 

0043 

POINTT=l 

0040 

POINTL=*0 

0045 

POINTH=l 

0046 

338 

CONTINUE 

0047 

DO  150  K000=i ,25 

0048 

DO  140  K0=0, 5 

0049 

READ(7'POINT2,END=77) (MAG (I) , I=l+K0*64 ,64+KO*64 ) 

0050 

140 

CONTINUE 

0051 

77 

CONTINUE 

0052 

DO  10  J=1 , 5 

005  3 

DO  5  1=1,127 

0054 

H2=H0 

0055 

H3=H1 

0056 

LCONS=I+l 

0057 

HO=AMMOUN (LCONS ,<J)  .AND. 2 55 

0058 

Hl=AMM0UN (LCONS ,  J+l) .AND. 255 

0059 

/  VERAG= (HO+H1+H2+H3) /4 

0060 

STOR ( I , 2 ) =0 

0061 

IF (AVERAG.LT. THRESH)  STOR(I,2)=l 

0062 

14 

CONTINUE 

0063 

5 

CONTINUE 

0064 

CALL  OUT (STOR, J+IJ) 

0065 

DO  756  IND=1, 128 

0066 

ST0R(IND,1)=ST0R(IND,2) 

0067 

7r'6 

CONTINUE 

006" 

10 

CONTINUE 

0060 

700 

CONTINUE 

0070 

POINT2=POINT2-l 

0071 

IJ=IJ+5 

0072 

150 

CONTINUE 

0073 

339 

FORMAT (32X, 15 (X, 14) ,/) 

0074 

341 

CONTINUE 

0075 

CALL  CLOSE (7) 

0076 

COORDI (l,6)=POINTH 

0077 

CALL  CHANGE 

0078 

CALL  REQUES (SEGMN) 

0079 

991 

CONTINUE 

0080 

END 
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When  the  Subroutine  out  is  called  it  uses  store  to  generate  Freeman-Code 
directions  for  a  row  and  three  such  rows  are  stored  in  the  array  number(128,3) . 

To  obtain  a  Freeman-Code,  the  patterns  of  Figure  13  are  assigned  the  indicated 
values  by  the  subroutine  out  and  stored  in  number(128,3) .  All  other  patterns 
are  assigned  a  -l. 


OOOl 

SUBROUTINE  OUT (STORE, J) 

0002 

INTEGER* 2  STORE(128,2) ,NUMBER(128,3) 

0003 

REAL* 4  ANGLES (15) 

0004 

REAL *4  THETA 

0005 

L0GICAL*1  FLAG 1, FLAG 2, FLAG 

0006 

DATA  ANGLES/-1. , 7. , 6. , 3. , 4. ,-l . ,5 . , 1 . ,-l. ,0 . , -1 . , 2 . , 3*-l ./ 

0007 

DATA  FLAG/. FALSE./ 

0008 

IF ( J . GT. 3) CALL  LOGIC (NUMBER, J) 

0009 

CALL  READEF ( 15 , LCODE) 

0010 

IF (LCODE . EQ . 2)  RETURN 

0011 

DO  10  1=1,127 

0012 

NUMBER (1,1) = NUMBER (1,2) 

0013 

NUMBER (1,2) = NUMBER (1,3) 

0014 

10 

CONTINUE 

0015 

DO  20  1=1,127 

0016 

INDEX=STORE ( 1 , 1 ) *2  *  *  3+STORE (1,2) *2  *  *2+ST0RE ( 1+1 , 1 ) *2  +  STORE (1+1,2) 

0017 

IF ( INDEX . EQ . 0 ) INDEX= 15 

0018 

NUMBER ( I , 3 ) = I INT (ANGLES ( INDEX) ) 

0019 

20 

CONTINUE 

0020 

RETURN 

0021 

END 

The  next  step  is  the  linking  of  directions  which  are  the  same,  and  appear 
sequentially  in  a  three-by-three  window.  Two  predicates  are  used  in  control¬ 
ling  the  statements  to  be  executed.  These  are 

p  =  The  element  of  number  is  a  continuation  of  a  segment  of  the  same 

direction. 

q  =  The  element  of  number  being  checked  is  continued  by  a  segment  of 
the  same  direction. 

The  cases  where  p=.true.  and  Q=.true.  are  illustrated  in  Figures  14a,  b. 

The  possible  predicates  and  the  actions  taken  when  the  predicates  are  true 
are  shown  in  Figure  14c. 

a.  paq 

Put  the  segment  number  of  the  line  that  the  element  of  number 
continues  into  first  15  bits  of  numberu^).  This  segement  number  is  extracted 
from  the  first  15  bits  that  numberu,2)  continues. 

b.  paq 

In  this  case  a  new  segment  needs  to  be  started.  The  starting  row 
and  columns  are  stored  in  the  array  lines,  pointl  contains  the  current 


segment  number  and  lines (pointl,1)  sets  the  row  number,  lines (po inti, 2) 
the  column  number  and  lines(pointl,3)  the  direction. 

c.  p  and  q 

This  happens  when  a  line  segment  of  the  same  direction  is  termin¬ 
ated.  The  action  taken  in  this  case  is  to  store  the  beginning  coordinates 
ending  coordinates,  and  direction  of  the  line  segment  in  coord i.  The  beginning 
point  of  the  segment  is  obtained  by  first  stripping  the  first  15  bits  off 
number ( i , 2 )  and  accessing  the  entry  of  lines  which  corresponds  to  this 
number.  This  gives  the  starting  point  of  the  segment,  while  the  final  point 
is  gotten  from  the  current  row  and  column  coordinates. 

d.  p  and  q 

It  is  an  isolated  direction  and  thus  it  is  stored  directly 

into  coordi. 


Subroutines  logic  and  guards  look  at  the  direction  numbers  in  number  and 
link  those  arrows  that  occur  sequentially  in  the  same  direction.  The 
arrays  used  to  do  the  bookeeping  at  this  stage  are  lines (,)  and  coordi (,) 
with  the  final  results  being  stored  in  coordi(,).  Long  polygonal  segments 
are  constructed  by  tracking  along  consistent  joins  of  these  line  segments 
at  each  point,  checking  for  possible  continuation  of  each  segment. 

The  subroutine  logic  drives  the  programs  which  produce  the  pseudo  Freeman 
code  and  do  the  bookeeping  functions.  The  data  from  which  it  computes 
line  segments  is  in  number(i28,3)  and  it  consists  of  the  Freeman  codes 
generated  from  the  last  three  rows  processed.  The  first  fifteen  bits  of 
number  are  used  to  store  the  segment  number  of  a  particular  entry. 

logic  processes  a  row  in  the  do  loops  of  line  fourteen  to  twenty-four. 

This  loop  begins  by  looking  to  see  if  the  element  number(i28,2)  is  a 
possible  edge  element,  and  if  It  finds  that  the  element  equals  "-1"  it 
looks  at  a  new  element  since  the  non-edge  elements  have  been  assigned  a  "-1." 
If  it  is  a  possible  edge  element  it  extracts  the  segment  number  and  calls 
guards  to  compute  the  values  of  p  and  q.  What  is  left  to  do  now  are  the 
actions  which  correspond  to  different  values  of  p  and  q  this  is  done  in 
lines  20  through  23. 

OOOl  SUBROUTINE  LOGIC (NUMBER, J) 

0002  INTEGER* 2  NUMBER(128, 3) 

0003  INTEGER* 2  LINES (500, 3) , COORDI (1500,6) 

0004  L0GICAL*1  P  J  IF  P=TRUE  THEN  ARROW  IS  A  CONTINUATION 

0005  LOGICAL*!  Q  I  IF  Q=TRUE  THEN  ARROW  IS  CONTINUED 

0006  INTEGER* 2  ROW,COL,POINTT,POINTH,POINTL, I 

0007  L0GICAL*1  INV , FLAG 

0008  INTEGER*2  SIGN,P0INSI 

0009  COMMON  INV, FLAG, SIGN 

0010  COMMON  LINES, POINTL,POINTH,POINTT, COL, 128 

0011  COMMON  /DTA/  COORDI 

0012  IF(J.EQ.l)  P0INTL=1 

0013  IF(J.EQ.l)  POINTH-1 

0014  DO  10  128=2,127 
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IF (NUMBER (128, 2) .EQ.-l) GOTO  10 
COL=NUMBER(I28,2) .AND. "177770 
CALL  GUARDS (NUMBER, P,Q, 128) 

ROW=NUMBER(I28,2) .AND. 7  !DIRECTION  CODE 

COL= (NUMBER (128, 2) .AND."177770)/8  ICOMPONENT  NUMBER 

IF(P.AND. .NOT.Q)  CALL  DUMP(J) 

IF ( ( .NOT.P) .AND.Q)  CALL  BEGIN (ROW, J) 

IF ( ( .NOT.P) .AND.Q)  NUMBER ( 128 , 2 )=NUMBER(I28, 2 ) .OR. (POINTL*8) 
IF (( .NOT.P) .AND. ( .NOT.Q) )  CALL  SINGLE (NUMBER, J, ROW) 

10  CONTINUE 
RETURN 
END 


logic  computes  the  values  of  p  and  q  by  using  the  subroutine  guards,  guards 
looks  in  a  three-by-three  neighborhood  of  states (128, 3)  and  computes  the 
value  of  p  and  q  for  states(i,2).  It  checks,  as  in  Figure  14,  for  the  appro¬ 
priate  values  of  p,  lines  12  -  57,  and  then  goes  on  to  compute  q,  line  60  -  end. 
This  processing  is  done  for  the  entire  file  with  the  final  results  being  stored 
into  coordi.  coordi (1500,6)  is  now  sorted  on  the  row  coordinates  of  its  ele¬ 
ments.  The  format  of  the  elements,  that  are  stored  in  coordi,  is  also  changed 
so  that  the  data  is  now 

CODE, BEG IN(j) .BEGIN (i) , DEL 


0001 

SUBROUTINE  GUARDS(STATES,P,Q,I) 

0002 

INTEGER*2  I, STATES (128,3) ,SEGMEU,CODE 

0003 

INTEGER* 2  CHOICP .CHOICQ 

0004 

LOGICAL* 1  P,Q,L 

0005 

P=. FALSE. 

0006 

Q=  .FALSE. 

0007 

L=. FALSE. 

0008 

CODE= (STATES (1,2). AND .7) 

0009 

INDEX=CODE+l 

0010 

SEGMEN- 0 

0011 

CHOICP=-l 

0012 

GOTO (10,20,30,40,50,60,70,80) INDEX 

0013 

10 

F( (STATES (1-1,2) .AND. 7) .EQ.O) CHOICP-STATES (1-1 , 2 ) 

0014  0014 

IF (CHOICP, EQ.-l)  GOTO  85 

0015 

SEGMEN=CHOICP . AND. "177770 

0016 

IF (SEGMEN, EQ.O) GOTO  85 

0017 

P= .  TRUE . 

0018 

STATES (1,2)— SEGMEN . OR . CODE 

0019 

GOTO  85 

0020 

20 

IF ((STATES (1+1,1) .AND. 7) . EQ . 1) SEGM£N= STATES (1+1,1) .AND. "177770 

0021 

IF (SEGMEN. EQ.O)  GOTO  85 

0022 

P= . TRUE . 

0023 

STATES (1,2) =  SEGMEN . OR . CODE 

0024 

GOTO  85 

0025 

30 

IF ( ( STATES 1,1) .AND. "177770. EQ. 2) SEGMEN-STATES (1,1). AND, "177770 

0026 

IF (SEGMEN. EQ.O)  GOTO  85 

0027 

P= . TRUE . 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 
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0 


P  and  Q 


P  and  Q 


P  and  0 


Figure  14.  Examples  of  geometric  relations  for  various 
predicates,  the  geometric  situation  when 
(a)  P  *  true,  and  (b)  Q  *  true. 
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STATES (1,2) =SEGMEN . OR . CODE 
GOTO  85 

40  IF ( (STATES (1-1,1) .AND. 7) . EQ. 3) SETMEN=  STATES (1-1,1) .AND. "177770 
IF(SEGMEN.EQ.O)  GOTO  85 
P= . TRUE . 

STATES (1,2) =SEGMEN . OR . CODE 
GOTO  85 

50  IF ( (STATES (1-1,2) .AND. 7) . EQ. 4) SEGMEN= STATES (1-1,2) .AND . "177770 
IF(SEGMEN.EQ.O)  GOTO  85 
P= . TRUE . 

STATES (1,2) =SEGMEN. OR. CODE 
GOTO  85 

60  IF ( (STATES (1+1,1) .AND. 7) . EQ . 5) SEGMEN=STATES ( 1+1 , 1) .AND. ”177770 

IF(SEGMEN.EQ.O)  GOTO  85 
P= . TRUE . 

STATES  (1,2)  =SEGMEN ..  OR.  CODE 
GOTO  85 

70  IF ( (STATES (1,1) . AND  .7) . EQ  .6)  CHO ICP= STATES ( 1 , 1) 

IF (CHOICP . EQ . -1 )  GOTO  85 
SEGMEN=CHOICP. AND. "177770 
P= . TRUE . 

STATES (1,2) =SEGMEN . OR . CODE 
GOTO  85 

80  IF( (STATES  (1-1,1) .AND. 7) .EQ. 7)  CHOICP=STATES (1-1,1) 

IF  (CHOICP . EQ .- 1)  GOTO  85 
SEGMEN=CHOICP. AND. "177770 
IF(SEGMEN.EQ.O)  GOTO  85 
P= . TRUE . 

STATES (1,2)  =  SEGMEN . OR . CODE 
GOTO  85 

85  CHOICP=-l 

CHOICQ=-l 

GOTO (100,200, 300 , 400 ,500,600, 700 , 800 ) INDEX 
100  IF( (STATES (1+1,2) .AND. 7) . EQ . 0 ) CHOI CQ=STATES (1+1,2) 

IF (CHOICQ . EQ. -1)  GOTO  850 
Q= . TRUE . 

GOTO  850 

200  IF ( (STATES (1-1,3) . AND. 7) .EQ. 1)  CHOICQ=STATES ( 1-1 , 3) 

IF (CHOICQ. EQ. -1  )  GOTO  850 
Q= . TRUE . 

GOTO  850 

300  IF ((STATES (I, 3) .AND. 7) .EQ. 2) CHOICQ=STATES (I , 3) 

IF (CHOICQ. EQ.-l)  GOTO  850 
Q= . TRUE . 

GOTO  850 

400  IF { (STATES (1+1,3) .AND. 7) . EQ . 3) CHOICQ=STATES (1+1 , 3) 

IF(CHOICQ. EQ.-l)  GOTO  850 
Q=. TRUE. 

GOTO  850 

500  IF ( (STATES (1+1,2) .AND. 7) .EQ. 4)  CHOICQ=STATES (1+1,2) 

IF (CHOICQ. EQ.-l)  GOTO  850 
Q= . TRUE . 
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0081  600  IF ( (STATES (1-1,3) .AND. 7) .EQ, 5)  CHOICQ= STATES (1-1,3) 

0082  IF (CHOICQ,EQ.-l)  GOTO  850 

0083  Q=.TRUE. 

0084  GOTO  850 

0085  700  IF ((STATES (I, 3) .AND. 7) . EQ. 6) CHOICQ=STATES ( I , 3) 

0086  IF (CHOICQ. EQ . -1)  GOTO  850 

0087  Q=.TRUE. 

0088  GOTO  850 

0089  IF ( (STATES (1+1,3). AND. 7) . EQ. 7)  CHOICQ=STATES ( 1  +  1 , 3) 

0090  IF (CHOICQ. EQ.-l)  GOTO  850 

0091  Q= . TRUE . 

0092  GOTO  850 

0093  850  RETURN 

0094  END 


SUBROUTINE  BEGIN (DIRECT ,J) 

INTEGER*2  LINES (500, 3) , COORDI (1500 , 6) 

INTEGER* 2  COL , I , J , DI RECT , POINTL , PO INTH , PO INTT 
LOGICAL* 1  INV, F  LAG 
INTEGER* 2  SIGN 
COMMON  INV, FLAG, SIGN 

COMMON  LINES , POINTL, POINTH , PO INTT, COL ,128 

COMMON  /DTA/COORDI 

POINTL=POINTL+l 

LINES (POINTL, 1)= 128 

LINES (POINTL, 2) =J 

L INES ( PO I NTL , 3 ) =  DI RECT 

RETURN 


0001  SUBROUTINE  DUMP(J) 

0002  INTEGER *2  L INES ( 500 , 3 ), COORDI ( 1500 , 6) 

0003  INTEGER* 2  I , J , DIRECT , POINTL, POINTH , POINTT 

0004  REAL* 4  BEGINX , BEGINY ,ENDX , ENDY 

0005  INTEGER*2  SIGN, COL 

0006  LOGICAL* 1  INV, FLAG 

0007  COMMON  INV , FLAG , S IGN 

0008  COMMON  LINES .POINTL , POINTH .POINTT, COL , 128 

0009  COMMON  /DTA/COORDI 

0010  DIRECT=LINES (COL, 3) 

0011  BEG INX*30+L INES (COL, 1 ) *7 

0012  BEG INY*  775-LINES (COL , 2 ) *  7 

0013  ENDX=30+( 128) *7 

0014  ENDY=775- (J) *7 

0015  GOTO (100, 200, 300, 400, 500, 600, 700, 800)  DIRECT* 1 

0016  100  ENDX-ENDX+7 

0017  GOTO  1000 

0018  200  BEGINX-BEGINX+7 

0019  BEGINY=BEGINY+7 

0020  GOTO  1000 

0021  300  BEGINY-BEGINY+7 

0022  GOTO  1000 
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0023 

400 

BEGINX-BEGINX-7 

0024 

BEGINY-BEG1NY+7 

0025 

GOTO  1000 

0026 

BEG I NX= BEG I NX -7 

0027 

GOTO  1000 

0028 

600 

ENDX=ENDX-7 

0029 

ENDY=ENDY-7 

0030 

GOTO  1000 

0031 

700 

EiiDY=ENDY-7 

0032 

GOTO  1000 

0033 

800 

ENDX-  ENDX + 7 

0034 

ENDY=ENLY  -7 

0035 

•GOTO  1000 

0036 

1000 

CONTINUE 

0037 

10 

COORDI (PCINTH ,  1)  *- (BEGINY-775 ) /7 

0038 

COCPOI (POINTH ,  2 )  = (BEGINX- 30 ) /7 

0039 

COORDI (POINTH , 3 ) =DIRECT 

0040 

20 

COORDI (POINTH ,4)=- (ENDY-775)/7 

0041 

COORDI  (P0TNTH,5)=  (E'.:DX-30)/7 

0042 

POINTH-PCINTH+1 

004  3 

RETURN 

0044 

END 

0001 

SUBROUTINE  CHANGE 

0002 

INTEGER*2  CoCRDi (1500,6)  !  INPUT  DATA 

0003 

INTEGER* 2  POINT  AMOUNT  OF  DATA 

0004 

LOGICAL*l  FLAG, BIO (2b) 

0005 

INTEGERS  DINES (500,3) ,POINTL, POINTH , POINTT 

0006 

INTEGER* 2  COL, I 28, SIGN 

0007 

DOGICAL*l  IN’; 

0008 

COMMON  INV , FLAG , SIGN 

0009 

COMMON  L INES , POINTL , POINTH , POINTT , COL . I 28 

0010 

COMMON  /DTA/COGRDI 

0011 

POINT-o 

0012 

FLA'i=  .  TRUE  . 

0013 

CALL  EP.RSET  (  37  ,  .TRUE.  ,  .FALSE.  ,  .FALSE  .  ,  .FALSE 

0014 

l')0 

CONTINUE 

C 

COORDI  FORMAT* CODE , TAIL ( J ) , TAIL (I) , DELTA 

0015 

DO  90  1=2, POINTH 

0016 

G’TO( 19, 20, 20, 20, 30, 40, 40, 40)  COORDI (1 ,3)+l 

0017 

10 

FLAG-COO  DT (1,2) .GE. COORDI (1,5) 

0018 

COORDI  (I  ,t>)  =  IABS  (COORDI  (I  ,2)  -COORDI  (1,5) 

0019 

IF ( . NOT . FLAG) COORDI (1,4) -COORDI (1,1) 

002'’) 

IF ( .NOT. FLAG) COORDI (I , 5! -COORDI (1,2) 

0021 

GOTO  86 

0022 

20 

FLAG-COORDI (I ,1) . GE . COORDI ( I , 4 ) 

0023 

COORDI (I ,6)-IABS (COORDI (I , 1 ) -COORDI ( I , 4 ) ) 

0024 

IF (FLAG)  COORDI (I ,4)=COORDI (I ,1) 

0025 

IF (FLAG)  COO PD T (I ,5) -COORDI (1,2) 

0026 

GOTO  88 

0027 

30 

FLAG=COORDI (1,5) . GE. COORDI ( I , 2 ) 

0028 

COORDI ( I , 6) = LABS (COORDI (1,2) -COORDI (1,5)) 

0029 

IF( .NOT. FLAG)  COORDI ( I ,4) -COORDI (I , 1) 

0030 

IF ( . NOT . FLAG )  COORDI (1,5) =COORDI (1,2) 

0031 

GOTO  88 

0032 

40 

FLAG=cOORDI (1,4) .GE. COORDI (1,1) 

0033 

COORDI (I , 6 ) = IABS (COORDI ( I ,1) -COORDI ( I , 4) ) 

0034 

IF (FLAG)  COORDI (1,4 ) -COORDI (1,1) 

0035 

IF (FLAG)  COORDI (1,5) =COORDI (1,2) 

0036 

GOTO  88 

0037 

88 

CONTINUE 

0038 

COORDI ( I , 1) =COORDI (1,3) 

0039 

COORDI (1,2) = COORDI (1,4) 

0040 

COORDI (I,3)=COORDI(I,5) 

0041 

COORDI (1,4) =COORDl (1,6) 

0042 

90 

CONTINUE 

0043 

CALL  SORT (POINTH , 2 ) 

0044 

DO  199  1=3 , POINTH- 3 

C 

WRITE (9,250)  COORDI (1,1) ,COORDI(I,2) ,COORDI(I,3) ,COORDI(I,4) 

004  5 

250 

FORMAT  (2X ,  11 ,  X ,  13 ,  X ,  13  ,  X ,  13  ) 

0046 

199 

CONTINUE 

0047 

RETURN 

0048 

END 

0001 

SUBROUTINE  SORT (POINT, KEY) 

0002 

INTEGERS  COORDI (1500 ,6)  !  TRANSMITTED  VIA  COMMON 

0003 

INTEGER*2  POINT  NUMBER  OF  ELEMENTS  TO  SORT 

0004 

INTEGER*2  KEY  )  WHICH  COLUMN  TO  SORT  ON 

000  5 

INTEGERS  L , R, K, X (4)  l  INTERMEDIATE  LOCATIONS 

0006 

INTEGER* 2  LINES (500, 3) 

000  7 

INTEGER*2  POINTL , POINTH , POINTT , COL , SIGN 

0008 

LOGICAL*l  FLAG.INV 

0009 

COMMON  INV, FLAG, SIGN 

0010 

COMMON  LINES, POINTL, POINTH, POINTT, COL, 128 

0011 

COMMON  /DTA/COORDI 

0012 

L=2 

0013 

R= POINT 

0014 

1 

CONTINUE  I  RETURN  POINT  FOR  OUTER  REPEAT 

0015 

DO  2  J=R , L , -1 

0016 

IF (COORDI (J, KEY) .GE. COORDI (J-l, KEY)  ) GOTO  2 

0017 

X(l)=COORDI (J-1,1) 

0018 

X (2 ) =COORDI (J-l ,2) 

0019 

X ( 3) =COORDI ( J-l , 3) 

0020 

X (4) = COORD I (J-l ,4 ) 

0021 

COORDI (J-1,1) -COORDI ( J , 1 ) 

0022 

COORDI ( J-l , 2 ) -COORDI ( J , 2 ) 

0023 

COORDI (J-l ,3) -COORDI (J,3) 

0024 

COORDI (J-l ,4) -COORDI (J,4) 

002  5 

COORDI (J ,1)=X(1) 

0026 

COORDI ( J , 2) =X (2 ) 

0027 

COORDI (J , 3 ) =X (3 ) 

0028 

COORDI (J,4)=X(4) 

0029 

K=  J 

0070 

2 

CONTINUE 

0031 

L=K+1 

0032 

DO  3  J=L,R 

0033 

IF (COORDI (J , KEY) . GE .COORD I (J-l , KEY) ) GOTO  3 

0034 

X (1) =COORDI ( J-l ,1) 

0035 

X (2) =COORDI (J-l, 2) 

0036 

X (3) =COORDI (J-l, 3) 

0037 

X (4) =COORDI (J-l ,4) 

0038 

COORDI ( J-l , 1 ) =COORDI (J,l) 

0039 

COORDI (J-l , 2) =COORDI (J,2) 

0040 

COORDI (J-l , 3) =COORDI (J,3) 

0041 

COORDI (J-l , 4) =COORDI ( J , 4) 

0042 

COORDI (J,1)=X(1) 

0043 

COORDI (J , 2 ) =X (2) 

0044 

COORDI (J, 3)=X(3) 

0045 

COORDI (J , 4) =X (4) 

0046 

K=J 

0047 

3 

CONTINUE 

0048 

R=K-1 

0049 

IF(L.LE.R)  GOTO  1 

0050 

END 

0001 

SUBROUTINE  SINGLE (STATES ,J , ROW) 

0002 

INTEGER* 2  I, J, DIRECT 

0003 

INTEGER*2  LINES (500 , 3) , COORDI (1500 ,6) 

0004 

LOGICAL*l  R,S,NR,NS,FLAG,  INV 

0005 

INTEGER*2  SIGN , ROW, COL, POINTT, POINTH,POINTL , DIR 

0006 

INTEGER* 2  CODE,  SEGMEN,CHOICP ,CHOICQ ,POINSI 

0007 

REAL *4  BEGINX , BEGINY , ENDX , ENDY 

0008 

INTEGER*2  STATES (128 , 3) 

0009 

COMMON  INV, FLAG, SIGN 

0010 

COMMON  LINES , PO INTL , POINTH , POINTT, COL ,128 

0011 

COMMON  /DTA/COORDI 

0012 

CODE= STATES (128,2) .AND. 7 

0013 

DIRECT=C-ODE 

0014 

BEGINX=30+ ( 128) *7 

0015 

BEGINY=775- (J) *7 

0016 

ENDX=30+ (138) *7 

0017 

ENDY=775- (J) *7 

0018 

GOTO (100,200,300,400, 500 , 600 , 700 , 800)  DIRECT+1 

0019 

100 

ENDX=ENDX+7 

0020 

GOTO  1000 

0021 

200 

BEGINX=BEGINX+7 

0022 

BEGINY=BEGINY+7 

0023 

GOTO  1000 

0024 

300 

BEGINY=BEGINY+7 

0025 

GOTO  1000 

0026 

400 

BEGINX=BEGINX-7 

0027 

BEGINY=BEGINY+7 

0028 

GOTO  1000 

0029 

500 

BEGINX=BEGINX-7 

0030 

GOTO  1000 

0031 

600 

ENDX=ENDX-7 

0032 

ENDY=ENDY-7 

0033 

GOTO  1000 

0034 

700 

ENDY=ENDY-7 

0035 

GOTO  1000 

0036 

800 

ENDX=ENDX+7 

0037 

ENDY=ENDY-7 

0038 

GOTO  1000 

0039 

1000 

CONTINUE 

0040 

COORDI (POINTH , 1)=- (BEGINY-775)/7 

0041 

COORDI (POINTH, 2)= (BEGINX-30) /7 

0042 

COORDI (POINTH , 3) =DIRECT 

0043 

COORDI (POINTH, 4 )=- (ENDY-775)/7 

0044 

COORDI (POINTH, 5)= (ENDX-30)/7 

0045 

COORDI (POINTH, 6 )=DIRECT 

0046 

P0INTH=P0INTH+1 

0047 

RETURN 

0048 

END 

The  fifth  coordinate  will  be  used  to  place  pointers  that  will  give  the  next 
piece  of  a  particular  polygonal  line,  if  there  is  one.  At  this  point  the 
scene  has  been  reduced  to  a  number  of  line  segments  of  different  length, 
each  having  one  of  eight  possible  directions.  The  next  step  is  to  link 
these  by  checking  to  see  if  there  is  a  possible  continuation  of  one  segement 
by  some  other  segment.  Such  a  linking  of  segments  is  the  function  of  the 
program  segarr.  To  begin  with,  all  elements  of  coordi (* ,5)  are  set  equal 
to  zero,  after  which  a  number  of  segments  is  built  up  in  the  following  steps: 

1.  Look  through  coordi(*,5),  and  if  an  entry  is  found  equal  to  zero 
then  proceed,  or  else  stop. 

2.  Start  a  new  segment  by  recording  the  location  of  the  zero  entry 
of  coordi (*,5)  in  SEGS. 

3.  Now  look  for  an  element  in  coordi  that  satisfies  linking  criteria  as 
given  in  Figure  16.  If  such  an  element  is  found,  two  different  cases  will  be 
considered.  Either  it  is  a  single  element,  or  it  is  a  segment  (more  than  one 
element).  The  two  alternative  courses  of  action  are: 

'  'Segment1  (a)  Link-up  data  structures,  as  in  Figure  15a,  which 
results  in  the  graphic  operations  (Figure  15c). 

'  'Single'  (b)  Link-up  data  structures,  as  in  Figure  15b,  which  re¬ 
sults  in  the  graphic  operations  (Figure  15d). 

After  these  segments  have  been  created  the  segment  list  is  looked  through, 
and  if  there  is  a  consistent  join  of  two  segments  whose  distance  apart  is 
less  than  three  units,  then  these  are  joined. 

The  program  segarr  links  together  the  segments  which  are  stored  in  the  common 
area  dta.  The  number  of  elements  in  coordi  is  passed  via  the  sixth  element 
of  coordi  (lines  9-ii).  coordi (*, 5)  will  be  used  to  store  the  pointers  and 
they  are  all  initialized  to  zero  in  lines  16  -  18.  The  line  with  label  85  is 
the  beginning  of  the  code  which  constructs  the  polygons  from  the  long  line 
segments.  First,  coordi(*,5)  Is  searched  for  a  zero,  i.e.,  a  segment  that 
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Figure  15.  Data  structures  (a),  (b)  and  operations  with  their  corresponding 
(c),  (d)  geometric  structures  and  geometric  operations. 
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hasn't  been  used  in  the  construction  of  a  polygon.  If  it  finds  a  zero*  it 
jumps  out  of  that  loop  and  begins  constructing  the  data  structure  which 
corresponds  to  a  polygon  (lines  20  -  23). 

The  pointer  which  keeps  track  of  the  numbers  of  polygons  points  gets  up¬ 
dated,  and  the  index  of  the  segment  to  be  processed  gets  stored  both  in 
segs (points ,i)  and  segs (points, 2) ;  also  an  "-1"  gets  stored  in  coordi(i,5) 
to  indicate  the  end  of  a  polygon  (lines  26  -  30).  The  call  to  where  com¬ 
putes  the  endpoint  of  the  current  segment  being  analyzed  and  stores  it  in 
the  array  search(3).  search(I)  containing  the  direction  of  the  segment 
and  search (2) ,  search(3)  the  column  and  row  coordinates  (lines  33  -  34). 

To  determine  if  there  is  a  segment  of  coordi  that  continues  the  segment 
at  pointv  a  call  to  find  is  made  on  line  34.  This  Subroutine  returns  an  o 
/  in  findp  if  there  is  no  continuation.  It  points  to  the  continuation  of  the 
segment  being  analyzed  if  one  has  been  found  that  passes  the  test  in  find 
(Figure  16,  Lines  19  -  86).  If  a  continuation  has  been  found  and  it  has 
coordi (findp, j)  =  o  then  it  is  one  of  the  original  long  line  segments  and 
it  is  added  to  the  list  being  constructed  (lines  35  -  40).  Alternatively 
it  may  be  that  there  is  a  continuation  of  the  element  being  tested  but 
that  this  continuation  is  a  segment.  In  this  case  the  list  which  is  being 
tested  is  added  to  the  continuation  list  (Lines  42  -  46)  with  the  corres¬ 
ponding  list  operations  as  in  Figure  15.  The  resulting  contour  has  many 
of  the  important  segments  in  it  but  there  are  many  gaps  in  the  contours 
(Figure  17)  which  should  be  closed.  One  obvious  method  of  closing  these 
is  to  search  through  all  the  segments  and  join  those  that  are  less  than 
a  certain  distance  apart  (Lines  52  -  77).  This  works  fairly  well,  as 


igure  18 

shows . 

0001 

INTEGER* 2  SEARCH(3)  1  CODE, HEAD  (J)  , HEAD ( I) 

0002 

INTEGER*2  COORDI (1500,6)  !  CODE , BEGIN (J) , BEGIN ( I) , DEL, POINTER 

0003 

INTEGER* 2  POINTH  !  ALIAS  FOR  POINT 

0004 

INTEGER* 2  SEGS (500,4)  )  BEGIN, END 

0005 

INTEGER* 2  POINTV , POINTS , PO INT , DI , DJ , F INDP 

0006 

INTEGER* 2  MIN  , INDEX, H,T, Y (3) 

0007 

LOGICAL  *1  BIO (25) ,FLAG, INTERN 

oooa 

COMMON  /DISOO/  SEGS, POINTS 

0009 

COMMON  /DT A/COORD I 

0010 

CALL  ERRSET (37, .TRUE. , .FALSE. , .FALSE. , .FALSE. ,31) 

0011 

POINT=COORDI (1,6) 

0012 

ISAVE=1 

0013 

K1=0 

0014 

602 

CONTINUE 

0015 

POINTS=0 

0016 

DO  90  1=2 , POINT+3  !  SET  ALL  POINTERS  TO  ZERO 

0017 

COORDI (I , 5) =0 

0018 

90 

CONTINUE 

0019 

85 

CONTINUE  !  RETURN  HERE  TO  BEGIN  A  SEGMENT 

0020 

DO  80  1= ISAVE+1 , POINT+1  !  SEARCH  FOR  UNUSED  ONES 

0021 

ISAVE= I 

0022 

IF (COORDI (1,5) .EQ.O)  GOTO  70 

0023 

80 

CONTINUE 

35 


TEST  1  -  IN -12 
TEST 2  -  JN-J2 
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TEST1 


Flo'jre  16,  'leometrlc  criteria  for  the  continuation  of 
a  polygonal  segment, 
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0024 

0025 

0026 

0027 

0023 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 

0061 

0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 


70  CONTINUE  !  JUMPS  OUT  OF  LOOP  HERE 
POINTV=I 

IF (POINTV. GT. (POINT) )  GOTO  20  !  FINISHED 

POINTS-POINTS+1  !  UPDATE  SEGMENT  POINTER 
SEGS (POINTS ,  1)  =  1  I  MARK  BEGINNING  OF  SEGMENT 
SEGS (POINTS, 2 )=I 
COORDI ( PO INTV , 5 ) = -1 

60  CONTINUE  1  RETURN  POINT  FOR  SEGMENT  CONSTRUCTION 
L=POINTV  1  INTERMEDIATE  STORAGE  LOCATION 

CALL  WHERE (COORDI (L,l) , COORDI (L, 2) ,COORDI(L,3) ,COORDI(L,4) .SEARCH) 
C  SEARCH  CONTAINS  ACTIVE  POINT  FOR  SEARCH 

CALL  FIND (SEARCH, FINDP , POINTV, NP)  !  FINDP  POINTS  TO  NEXT  OR  0 

IF (FINDP . EQ .0)  GOTO  30  !  CANNOT  CONTINUE 

IF (COORDI (FINDP, 5) .NE.O)  GOTO  50 

COORDI (POINTV , 5) =FINDP 

SEGS (POINTS , 2 ) -FINDP 

COORDI (FINDP, 5 )=-l 

PO INTV=F INDP 

GOTO  60 

50  CONTINUE!  MERGE  LISTS 

COORDI (POINTV , 5 ) -SEGS (NP , 1 ) 

SEGS (NP , 1 ) =  SETS ( POINTS , 1 ) 

POINTS- POINTS-1 

GOTO  BS  !  BEGIN  A  NEW  SEGMENT 

COORDI (POINTV ,  5 )  — 1 

IF ( SEGS ( PO INTS , 1 ) . EQ . SEGS ( PO INTS , 2 ) ) COORD I ( POINTV , 5 ) =  0 
IF (SEGS (POINTS, 1) .EQ. SEGS (POINTS, 2) )  POINTS=POINTS-l 
GOTO  85 

20  CONTINUE 

K1=K1+1 
MIN=5000 
H=SEGS (K1 ,2) 

IF(H.EQ.-l)  GOTO  110 

CALL  WHERE ( COORDI (H,l) ,COORDI(H,2) ,COORDI(H,3) .COORDI (H, 4) ,Y) 

DO  150  J=l, POINTS 
T=SEGS (J ,1) 

IF (T.EQ.-l)  GOTO  150 
DJ= IABS (COORDI (T, 3) -Y (3) ) 

DI= IABS (COORDI (T , 2 ) -Y ( 2 ) ) 

IDIS=DJ+DI 

IF (ID IS . LT.MIN) INDEX=J 
IF(IDIS.LT.MIN)  MIN-IDIS 
150  CONTINUE 

IF (MIN. GT. 3) GOTO  110 
ICONST-SEGS ( INDEX , 2 ) 

IF (SEGS ( INDEX, 1) .EQ. COORDI (ICONST, 5) )  GOTO  110 
COORDI (SEGS (Kl, 2) , 5) -SEGS (INDEX, 1) 

SEGS ( K1 , 2 ) -SEGS ( INDEX , 2 ) 

IF (Kl.EQ. INDEX)  GOTO  110 
SEGS  ( INDEX,  1)— 1 
SEGS  (INDEX, 2 )— 1 
KL-Kl-1 

110  CONTINUE 
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0076 

IF (K1 .NE .POINTS) GOTO20 

0077 

21 

CONTINUE 

0078 

CALL  INITT(160) 

0079 

CALL  DWINDO(-50. ,1000. ,10. ,850.) 

0080 

CALL  CHRSIZ (3) 

0081 

CALL  COMPLE(LX) 

0082 

CALL  DISPLA 

0083 

CALL  CHRSIZ (3) 

0084 

CALL  FINITT(0, 780) 

TYPE  *,'$LX',LX 

C 

ACCEPT* , IJO 

c 

IF (IJO ,EQ .0)  STOP 

0085 

230 

FORMAT ( 3X , 6 ( 15 , X) ) 

0086 

LTIME=1 

0087 

670 

CONTINUE 

0088 

TYPE  *,'*?' 

0089 

READ (5 ,222) IJO 

0090 

IF (IJO,EQ.O)  STOP 

0091 

IF (LTIME . EQ . 1) CALL  INTTT(160) 

0092 

LTIME=TIME+1 

0093 

CALL  DWINDO(-50. ,1000. ,10. ,850.) 

0094 

CALL  CHRSIZ (3) 

0095 

CALL  CHRSIZ (3) 

0096 

CALL  DISPL1 ( IJO) 

0097 

CALL  CHRSIZ (3) 

0098 

CALL  F1NITT( 0,780) 

0099 

IF (IJO .NE. 0)  GOTO  670 

0100 

222 

FORMAT (13) 

0101 

CALL  CHRSIZ (3) 

ccccc 

DO  500  1=1, POINTS 

C  IF  (SEGS  ( 1 , 4)  .  LT  .20)  GOTO  500 

C  IF (SEGS (1,1) .EQ.-l) GOTO  500 

C  IHEAI>  SEGS  (1,1) 

•  C  L=SEGS (1,2) 

C  CALL  WHERE (COORDI(L,l) ,C00RDI(L,2) ,COORDI(L,3) ,COORDI(L,4) , SEARCH) 

C  ICYCLE=TABS (COORD I ( IHEAD , 3 ) -SEARCH ( 3 ) ) 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 


SUBROUTINE  FIND (ACTIVE ,FINDP ,POINTV,NP) 

INTEGER* 2  ACTIVE (3)  l  POINT  FROM  WHERE  SEARCH  IS  MADE 
INTEGER* 2  POINT 
INTEGER* 2  TEST1 ,TEST2 

INTEGER* 2  FINDP  l  INDEX  OF  POINT  FOUND  OR  ZERO 
INTEGER* 2  COORDI (1500 ,6)  !  DATA  TO  BE  SEARCHED 

INTEGER* 2  SEGS (500,4)  !  SEGMEN  POINTERS 

INTEGER* 2  DI,DJ, DEL, PO INTV 

INTEG£R*2  NP  .*  INDEX  OF  SEGS  FOR  HEAD  OF  MERGE 

INTEGER*2  POINTS 

COMMON  /DTA/COORDI 

COMMON  /DISOO/SEGS, POINTS 

POINT=COORDI (1,6) 

FINDP- 0 
JDIS-1 
JLIM-0 
J2-ACTIVE (2) 
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I2=ACTIVE(3) 

DO  5  L=100,l,-1 
IF ( (POINTV+L) . LT.l)  GOTO  4 
IF (POINTV+L.GT. POINT) GOTO  4 
JN=COORDI ( PO INTV+L , 2 ) 

IN=COORDI (POINTV+L, 3) 

TEST1= IN-12 

IF (IABS (TESTl) . GT . JDIS) GOTO  4 
TEST2=JN-J2 

IF (IABS (TEST2) .GT. JDIS)  GOTO  4 
GOTO (10,20,30,40,50,60,70,80)  ACTIVE (1) +1 
10  IF (TESTl, LT.JLIM) GOTO  4 

GOTO  90 

20  IF((TEST1, LT.JLIM) .AND, (TEST2.GT.JLIM))  GOTO  4 

GOTO  90 

30  IF (TEST2 . GT . JLIM)  GOTO  4 

GOTO  90 

40  IF ( (TESTl. GT. JLIM) .AND. (TEST2.GT. JLIM))  GOTO  4 

GOTO  90 

50  IF (TESTl. GT. JLIM)  GOTO  4 

GOTO  90 

60  IF ( (TESTl. GT. JLIM) .AND. (TEST2. LT.JLIM))  GOTO  4 

GOTO  90 

70  IF (TEST2.LT. JLIM)  GOTO  4 

GOTO  90 

80  IF ( (TESTl. LT.JLIM) .AND. (TEST2.LT. JLIM))  GOTO  4 
GOTO  90 

90  CONTINUE 

FINDP=POINTV+L 

IF (COGRDI (POINTV+L+5) .EQ.O)  RETURN  l  A  SIMPLE  CONSTRUCT 
DO  1  J=1 , POINTS-1  I  IS  CANDIDATE  THE  HEAD  OF  A  LIST 
NP=J 

IF  SEGS(J,1) .EQ.FINDP.AND. (COORDI (SEGS(J,2) ,5) .EQ.-l) )  RETURN 
1  CONTINUE 

4  FINDP=0 

IF ( (POINTV-L) ,LT. 1)  GOTO  5 

IF ( (POINTV- L-l. ) .GT. POINT)  GOTO  5 

JN=COORDI (POINTV-L , 2 ) 

IN=COORDI (POINTV-L, 3) 

TEST=IN-12 

IF (IABS (TESTl) .GT. JDIS)  GOTO  5 
TEST2=JN-J2 

IF (IABS (TEST2) .GT. (JDIS) GOTO  5 

GOTO (100,200,300,400,500,600,700,800)  ACTIVE ( 1 ) +1 
100  IF  (TESTl .LT.JLIM)GOTO  5 

GOTO  900 

200  IF (TESTl. LT.JLIM) .AND. (TEST2.GT. JLIM)  GOTO  5 

GOTO  900 

300  IF (TEST2 . GT. JLIM)  GOTO  5 

GOTO  900 

400  IF((TEST1. GT. JLIM) .AND. (TEST2.GT. JLIM))  GOTO  5 

GOTO  900 

500  IF (TESTl. GT. JLIM)  GOTO  5 
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! 

i 

J 


0071 

GOTO  900 

* 

0072 

600 

IF ( (TESTl . GT . JLIM) .AND. (TEST2 . LT. JLIM) ) 

GOTO  5 

* 

0073 

GOTO  900 

'' 

0074 

700 

IF (TEST2.LT. JLIM)  GOTO  5 

0075 

GOTO  900 

0076 

800 

IF ( (TESTl . LT .JLIM) .AND . (TEST2 . LT .JLIM) ) 

GOTO  5 

0077 

GOTO  900 

i 

0078 

900 

CONTINUE 

0079 

FINDP=POINTV-L 

0080 

IF (COORDI (POINTV-L, 5) . EQ . 0)  RETURN  !  A 

SIMPLE  CONSTRUCT 

0081 

DO  2  J=1 ,P0INTS-1 

0082 

NP=J 

0083 

IF (SEGS ( J , 1 ) . EC- FINDP. AND. (COORDI (SEGS(J 

,2) ,5) .EQ.-l)) RETURN 

0084 

2 

CONTINUE 

0085 

FINDP=0 

0086 

5 

CONTINUE 

0087 

RETURN 

0088 

END 

This  process  results.  In  most  cases.  In  a  closed  curve  that  can  be  analyzed 
by  the  Fourier  Classification  process  but  there  are  also  cases  where  the 
contours  produced  are  not  suitable  for  processing  but  must  be  first  modified 
by  an  operator  (Figure  4,  5)  before  they  can  be  used.  In  this  case,  the 
next  step  In  the  processing  allows  an  operator  to  Interactively  modify  the 
contours  produced  so  that  they  are  closed  curves  and  can  be  analyzed  by  the 
Fourier  descriptor  programs.  Programs  that  are  used  to  modify  the  computer¬ 
generated  polygons  are  documented  In  Reference  2.  Examples  of  how  they  work 
are  in  Figures  3,  4,  19  and  20. 

FOURIER  DESCRIPTOR  METHODS 


The  library  that  is  stored  in  the  computer  does  not  use  the  (x,y)  coordinates 
of  the  polygons,  but  first  transforms  them  via  the  Fast  Fourier  Transforms 
and  stores  the  "Fourier  Descriptors."  These  Fourier  Descriptors  are  defined 
as  follows:  A  closed  curve  can  be  thought  of  as  a  function  of  a  complex  var¬ 
iable,  z(t),  parametrized  by  arc-length  t.  We  can  normalize  and  have  the 
curve  described  by  z(t),  o  <  t<  2*pi.  If  we  go  around  the  contour  more  than 
once,  we  get  a  periodic  function,  which  can  be  expanded  in  a  convergent  Four¬ 
ier  series.  The  Fourier  Descriptor  of  the  curve  is  defined  to  be  the  Complex 
Fourier  series  expansion  of  z (t)  which  is  given  by  the  formula 
00 

z(t)  =  l  A(n)  eint  where 

A<«*>  -  0/2U  Z(t)  e“intdt 

(See  Figure  21,  22) 

Thus  the  Fourier  Descriptors  (32  of  the  a (i) )  for  each  element  in  the  library 
are  computed  and  stored  Into  memory.  The  contour  of  the  unknown  plane  is  then 
found,  the  Fourier  coefficients  for  this  unknown  are  calculated,  and  the  angu¬ 
lar  data  necessary  Is  obtained  by  finding  the  element  i  the  library  whose 
Fourier  coefficients  are  closest  to  the  unknowns. 
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tlons  entered  by  the  operator. 


Figure  20.  New  polygons  obtained  from  corrected  polygons 
of  Figure  19. 
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Flnure  22.  Three  basic  neonetrlcal  shades  and  their 
Fourier  coefficients. 
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In  order  to  compute  the  Fourier  Transforms,  a  closed  curve  description  of 
the  target  to  be  analyzed  has  been  produced  by  the  computer  and  the  operator. 
The  description  of  this  curve  consists  of  a  sequence  of  x,y  coordinates,  which 
are  the  vertices  of  a  polygon.  As  the  first  step  the  length  of  this  contour 
is  computed,  and  the  contour  is  resampled  at  a  spacing  chosen  to  make  the  total 
number  of  samples  a  power  of  two.  This  polygon  is  then  filtered  to  remove 
noise,  and  the  Fourier  descriptor  is  computed  by  taking  the  Fast  Fourier  Trans¬ 
form  of  this  sequence  of  (x,y)  coordinates. 

If  two  polygons  are  congruent  in  the  plane  then  they  can  be  shown  to  be  so  by 
a  sequence  of  rotations,  translations  and  contractions  followed  by  a  point-to- 
point  comparison.  If  we  have  two  congruent  triangles  represented  by  a  sequence 
of  x,y  coordinates  they  can  be  shown  to  be  congruent  by  first  rotating  both  so 
that  their  longest  side  lies  on  the  x  axis,  doing  separate  contractions  so  that 
they  both  have  the  same  area  then  doing  a  point-by-point  comparison  starting  at 
the  greatest  x  coordinate.  It  is  clear  that  the  point-by-point  comparison  must 
be  done  starting  at  the  same  place  on  both  triangles,  and  continuing  at  equi- 
distantly  sampled  points  in  order  for  this  process  to  be  meaningful.  The  geo¬ 
metric  transformations  used  to  show  two  polygons  congruent  translate  into  the 
frequency  domain  as  shown  in  Table  I. 


TABLE  I.  EQUIVALENT  OPERATIONS 


TIME  DOMAIN 
Translation 
Rotation 

Comparison  point 
change 


FREQUENCY  DOMAIN 
Addition  to  a(0) 

Multiplication  of  series  by  a  constant 
Multiplication  of  a(j)  by 
exp(ijt) 


In  order  for  a  comparison  to  be  meaningful  in  the  frequency  domain,  a  "normali¬ 
zation"  in  the  frequency  domain  must  be  done  similar  to  the  geometric  normali¬ 
zation  that  has  been  done  for  the  triangles.  This  normalization1  must  be  done, 
using  only  the  operations  which  are  listed  on  the  right  of  Table  I. 

First,  a (O)  is  set  equal  to  zero  to  normalize  position.  Size  normalization  is 
accomplished  by  dividing  each  coefficient  by  the  absolute  value  of  ad).  To 
normalize  the  point  where  the  comparison  is  to  begin,  we  require  that  the  phase 
of  the  two  coefficients  of  largest  magnitude  be  zero.  For  a  simple  closed 
curve  that  does  not  cut  itself  ad)  is  the  coefficient  of  largest  magnitude. 
Some  polygons  and  their  normalized  Fourier  coefficients  appear  in  Figure  22. 

Let  a (L)  and  a(K)  be  two  non-zero  coefficients  of  the  Fourier  series.  The 
normalization  multiplicity  of  the  coefficients  a(L)  and  a(io  is  defined  to  be 
M=abs (k-l) .  Some  of  the  geometric  significance  of  m  in  the  case  where  ad)  and 
a(L)  are  the  only  non-zero  coefficients  is  given  by  the  following  proposition: 

Let  z(t)  =  Ad)*expdt)  +  A(L)  *exp  (iLt)  with  abs(Ad))  >  abs(A(L))  >  0 


"d  P*  Wallace  and  0.  R.  Mitchell,  "Local  and  Global  Shape  Descrip¬ 
tion  of  Two  and  Three  Dimensional  Objects,"  School  of  Electrical  Engineering, 
Purdue  University,  September  1979. 
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PROPOSITION  1.  If  abs (A(i) )  >  abs (l*a (l) )  then  the  closed  curve  described 
by  z(t)  has  no  intersections. 

PROPOSITION  2.  If  abs (A(i) )  =  abs (l*a (l) )  then  the  closed  curve  described 
by  z(t)  has  m  =  absd  -  l)  cusp  points.  The  angles  at  these  points  are 
convex  if  l  <  o  and  concave  if  l  >  o  (see  Figure  23). 

PROPOSITION  3.  The  function  abs(z(t))  has  m  maximum  points  and  m  minimum 
points. 

Let  z ( t)  =  \  A( j)exp (ijt) 

j  =  —00 

PROPOSITION  4.  The  requirement  that  a(L)  and  a(K)  have  zero-phase  angle 
can  be  satisfied  by  m  different  orientation/starting  point  combinations. 

PROPOSITION  5.  m*x  | a ( j ) | ^  Z  ™oo  =  ! a (l) |  if  the  associated  curve  has  no 

intersections.2 

Thus  for  a  figure  whose  second  greatest  coefficient  is  a(L),  there  are 
m  =  abs (L  -  l)  possible  ways  to  normalize  this  figure.  In  order  for  an 
accurate  comparison  to  be  possible  the  normalization  chosen  for  like  figures 
wust  be  the  same.  We  use  the  following  method  to  choose  the  normalization: 

1.  Calculate  the  Fourier  coefficients  for  the  m  possible  normalizations. 

2.  For  each  of  the  m  normalizations  calculate 

y  re(a.)  la. ! 

L  x  x 1 

3.  Use  the  normalization  which  maximizes  the  above  quantity. 

The  pattern  recognition  method  begins  by  constructing  a  three-dimensional 
representation  of  the  target  to  be  analyzed.  A  library  of  polygons  is  then 
constructed  which  are  the  projections  of  the  three-dimensional  object,  as 
seen  from  different  views.  From  this  library  of  polygons  a  library  of 
Fourier  Descriptors  Is  computed  using  normalization  described  above,  and 
stored  Into  the  computer.  When  an  unknown  is  to  be  analyzed  and  its  contour 
is  found,  the  Normalized  Fourier  Descriptors  are  calculated  and  these  num¬ 
bers  are  compared  to  the  library  entries  via  the  difference. 

I  I ai  ~  LIB.(i) | 


2.  K.  Phillips  and  R.  Machuca,  "The  Geometry  of  Closed  Curves  Parame- 
tized  by  Fourier  Series,"  Research  Memorandum,  White  Sands  Missile  Range, 
Instrumentation  Directorate,  Advanced  Technology  Office. 


Once  the  closest  element  of  the  library  Is  found,  an  Interpolation3  Is  done 
to  get  an  accurate  measure  of  the  aspect  angles.  When  this  is  done  the  pro¬ 
cessing  of  the  frame  Is  finished. 

We  have  described  an  interactive  system  that  could  be  used  to  obtain  aspect 
angle  information  from  one  frame  of  video.  Before  this  system  can  be  made 
completely  automatic  there  must  be  research  done  with  regards  to  two  diffi¬ 
cult  problems.  One  is  the  automatic  choosing  of  a  threshold  which  would 
separate  possible  target  points  from  the  background.  An  approach  being  con¬ 
sidered  is  an  adaptive  procedure  for  choosing  p  where  p  would  be  incremented 
if  the  size  of  the  ellipse  defined  by  the  coefficients  a(l)  and  a(-l)  increased; 
and  p  would  be  decreased  if  the  area  of  this  ellipse  decreased.  The  other  is 
the  extraction  of  the  target  from  the  polygonal  representation  of  the  scene. 

Both  are  difficult  problems  which  will  require  much  research  before  a  satis¬ 
factory  solution  can  be  found. 


3.  T.  P.  Wallace  and  0.  R.  Mitchel,  ibid. 
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Appendix  A 

FINDING  EDGES  IN  NOISY  SCENES 


Research  into  methods  of  Identifying  edges  in  a  noisy  scene  has  been  an 
active  field  of  investigation  for  many  years.  Treatment  of  the  subject  may 
be  found  in  many  books  written  over  the  past  decade4'5'6  and  many  different 
approaches  are  proposed.  Recently  a  survey  and  comparative  analysis  of  the 
methods  was  made.7 

The  body  of  this  appendix  is  segmented  into  four  parts.  In  the  first,  we 
derive  and  define  a  "Moment  Operator"  which  we  show  to  work  well  for  step 
and  ramp  edges.  Then,  we  define  and  characterize  second  order  edges  using 
the  concept  of  the  rotation  of  a  point  in  a  vector  field  and  develop  the 
detector  analytically.  In  Section  3  we  develop  the  algorithms  for  imple¬ 
menting  the  previously  defined  operators.  Finally,  in  Section  4,  these 
algorithms  are  evaluated  using  ROC  curves  and  compared  with  previously 
known  techniques. 

The  detection  of  edges  to  isolate  objects  in  a  scene  is  motivated  by  many 
distinct  problems.  One  such  problem  arises  in  a  tracking  system  where 
the  input  video  Image- is  analyzed  and  the  object  to  be  tracked  identified. 
Subsequent  input  and  feedback  to  the  drive  controls  causes  the  sensor  to 
re-orient  to  a  new  position  in  an  attempt  to  maintain  the  same  x-y  coordi¬ 
nate  position  for  the  object  in  the  field  of  view.  While  this  problem 
motivated  the  research  that  led  to  this  paper,  the  results  herein  discussed 
are  much  broader  in  scope  and  application.  The  constraints  imposed  by  this 
problem  led  to  a  method  that  is  useful  in  high  data  throughout  systems. 

Section  1.  Edges  from  Moments 

First  order  edge  detection  methods  work  in  the  following  way:  A  picture 
function  f(x,y)  is  transformed  to  another  picture  function  F(x,y)  =  Tf(x,y) 
in  such  a  way  that  the  edges  of  objects  in  the  scene  will  be  in  the  set 
{  (x,y):F(X,y)  >  W}  for  some  W.  The  usual  method  is  to  transform  the  picture 
using  T  equal  To  the  gradient  operator.  Different  edge  detection  methods 
correspond  to  different  numerical  approximations  to  the  gradient. 


4.  A.  Rosenfeld  and  A.  Kak,  "Digital  Picture  Processing,"  Academic 
Press,  New  York,  New  York,  1976 

5.  B.  Lipkin  and  A.  Rosenfeld,  "Picture  Processing  and  Psychopictorics" , 
Academic  Press,  New  York,  New  YOrk,  1970 

6.  W.  Pratt,  "Digital  Image  Processing,"  John  Wiley  and  Sons,  New  York, 
New  York,  1978 

7.  I.  Abdou,  "Quantitative  Methods  of  Edge  Detection,"  Image  Processing 
Institute,  University  of  Southern  California,  Los  Angeles,  California,  1978 
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The  method  used  in  our  edge  detection  program  is  not  based  on  derivatives. 
To  reduce  the  effect  of  noise,  this  edge  detection  method  uses  integrals. 

The  reasoning  for  the  use  of  moments  to  find  edges  is  as  follows: 

■  A  digitized  picture  can  be  thought  as  a  lamina  whose  density  at 
each  point  is  f(x,y),  so  points  of  high  intensity  correspond  to  points  of 
high  density. 

‘  A  point  (a,b)  on  an  edge  in  the  original  function  (see  Figure  A-l) 
would  correspond  to  a  point  in  this  lamina  (digitized  picture)  with  high 
densities  on  one  side  and  lower  densities  on  the  other  side. 

Thus  if  we  look  at  a  small  lamina  centered  at  point  (a,b)  and  compute  the 
center  of  mass  of  this  small  lamina,  we  can  expect  the  center  of  mass  to 
lie  within  an  area  of  high  densities. 


3TTTTr('ii 

jj||  j,:  =  regions  of  high  density 


Figure  A-l.  Example  center  of  mass  vectors  for  (1)  and  edge 
and  (2)  a  region  of  uniform  intensity. 


Suppose  we  now  look  at  a  point  (c,d)  such  that  the  densities  around  it  are 
fairly  constant.  Then  the  center  of  mass  of  a  small  lamina  about  it  would 
be  close  to  (c,d).  In  this  case,  a  vector  from  (c,d)  to  the  center  of  mass 
would  be  very  small  compared  to  a  vector  from  (a,b)  to  the  center  of  mass 
In  the  previous  case. 

We  conclude  that  one  way  to  transform  f(x,y)  to  F(x,y)  such  that  edges  of 
the  original  picture  lie  in  the  set  F(x,y)  >  W  is  to  replace  every  f(x,y) 
by  the  length  of  the  vector  from  (x,y)  to  tTFe  center  of  mass  of  a  small 
lamina  centered  about  (x,y).  That  is,  F(x,y)  is  the  magnitude  of  the  vector 
from  (x,y)  to  the  center  of  gravity  of  a  square  lamina  centered  at  (x,y) 
whose  density  Is  given  by  the  picture  function  f(x,y). 

Figure  A-2(b)  is  an  example  of  how  this  method  works  on  a  scene  (Figure  A-2(a)) 
typical  of  those  we  study  at  WSMR. 


[a)  IMAGE  OF  ROCKET  AND  PLUME.  THE 
PLUME  IS  THE  LARGE  REGION  OF 
HIGHEST  INTENSITY. 


(c)  THE  VECTOR  FIELD  GENERATED  BY 
THE  MOMENT  OPERATOR. 


(b)  RAMP  AND  STEP  EDGES  FOUND  BY 
USING  THE  MOMENT  OPERATOR. 
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(d)  SECOND  ORDER  EDGES  DETECTED 
BY  USING  THE  VECTOR  FIELD. 


Figure  A-2.  Rocket  and  results  of  processing. 


Once  the  coordinates  (Y,Y)  of  the  center  of  mass  of  a  lamina  about  (x,y) 
are  calculated,  the  direction  of  the  edge  (if  any)  can  easily  be  found. 
Since  (7,7)  points  to  where  the  intensity  of  the  picture  is  the  highest, 
the  direction  of  the_edge  is  perpendicular  to  the  direction  of  the  vec¬ 
tor  from  (x,y)  to  7,Y.  If_we  take  (x,y)  =  (0,0),  then  the  direction  of 
the  edge  is  e  =  Arctan  (Y./X)  +  ~72. 

Thus  this  model  gives  for  each  point  in  the  scene  a  quantity  that  meas¬ 
ures  the  probability  that  a  point  is  an  edge  point  and  a  direction  which 
is  the  direction  of  a  possible  edge  through  that  point. 

The  model  introduced  in  Section  I  will  not  work  for  roof  edges.  This  is 
because  at  the  very  peak  of  the  roof,  exactly  where  the  edge  is  situated, 
both  X  and  Y  are  equal  to  zero.  In  order  to  detect  roof  edges  we  need 
to  take  advantage  of  the  direction  information,  and  as  Figures  6(a),  (b) 
and  (c)  show  we  need  to  detect  the  shearing  cause  by  the  change  in  direc¬ 
tion  of  the  vector  field  at  the  edge  points.  One  way  of  doing  this  is  by 
using  a  tool  from  the  theory  of  vector  fields,  namely  the  rotation  of  a 
vector  field  about  a  point. 


If  a  curve  r  on  the  plane  (scene)  is  given  in  the  form 
r:  x  =  x(t),  y  =  y(t)  a  <  t  <_  b 

then  0(t)  =  { <J>[x( t) ,  y(t)],  ijj[x(t),  y(t)]  is  defined  on  the  interval  [a,b] 
(see  Figure  A-3). 


Figure  A-3.  A  curve  r  and  its  corresponding  vector  field  o(t) 


For  each  t  e[a,b]  there  is  determined  an  angle,  the  angle  in  radians  between 
O(t)  and  0(a)  measured  from  0(a)  to  0(t).  This  angle  is  a  many-valued  func¬ 
tion  (vanishing  for  t  *  a)  Is  designated  by  e(t)  and  called  an  angular 
function  of  the  field  o  on  a  curve  r.  The  rotation  of  the  field  0  on  the 
curve  r  is  defined  to  be 

Y(*,r)  =  ^  [0(b)  -  ©(a)] 
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If  r  is  a  closed  Jordan  curve,  then  the  rotation  is  found  by  subdividing 
r  into  two  curves  (not  closed),  computing  the  rotation  of  each,  and  adding. 
In  the  following,  r  is  taken  to  be  a  small  circle  about  a  point. 

We  can  write  the  rotation  as 

V  =  ^  [0(b)  -  «.)]  =  k  it. 

With  0(t)  =  Arctan  Y/Y  +  -rr/2,  we  make  the  following  observations: 

If  0(t)  =  constant,  then  =  0  and  y  =  0.  So  y  =  0  when 

x  =  a  point  on  the  edge  of  an  object  in  a  scene  (see  Figure  A-4). 


Figure  A-4.  Vector  Field  at  a  step  or  ramp  edge  point. 


Section  2.  Second  Order  Edges 

After  a  scene  is  processed  by  the  moment  edge  detector,  each  point  is 
assigned  a  direction  and  a  magnitude.  In  effect  this  specifies  a  vector 
at  each  point  of  the  plan  In  question;  l.e.,  these  vectors  define  a  vector 
field  over  the  scene.  An  Important  tool  in  the  study  of  vector  fields  is 
the  rotation  of  a  vector  field.  8,9  To  define  the  rotation  of  a  vector  field, 
suppose  a  vector  of  the  vector  field  $  at  the  point  (x,y)  Is  given  by 

®(x,y)  =  U(x,y),  ^(x,y)} 

4>(x,y)  =  I(x,y) 

t|>(x,y)  =  T(x,y) 

If  ©  is  symmetric  about  x  and  r  is  a  small  circle  about  x  =  edge  point 
on  a  roof  edge,  see  Figure  5,  then  write  r  *  r,  +  r?  (where  r.,  -  one  half  of 
the  circle  and  =  the  other  half)  . 


8.  J.  Milnor,  "Topology  from  the  differentiable  viewpoint,"  University 
Press  of  Virginia,  Charlottesville,  Virginia,  1965 

9.  A.  H.  Stroud,  "Approximate  Calculation  of  Multiple  Integrals," 
Prentice  Hall,  Englewood  Hills,  New  Jersey,  1971 
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Figure  A-5.  Vector  field  at  a  roof  edge  point. 


Figure  A-6(a)  and  A-6(b)  are  examples  of  how  these  observations  can  be  used 
to  detect  second  order  edges. 


Section  3.  Algorithms  for  Implementation 

a.  Calculation  of  Moment.  Since  we  are  interestedjn  real  time  appli¬ 
cations  of  these  methods  we  simplify  the  calculation  of  X  and  Y  by  setting 

M  =  J^h  J^k  f(x  +  t,  y  +  u)  dt  du  =  1 

This  can  be  justified  by  observing  that  M/4hk  is  the  average  of  the  inten¬ 
sities  over  a  small  neighborhood  of  (x,y)  and  so  this  value  can  be  approxi¬ 
mated  by  the  average  value  of  intensities  over  the  entire  picture.  This 
would  then  be  just  a  scale  factor  and  so  could  be  left  out. 

To  calculate  the  integrals  involved  we  use  an  integral  formula10  of  order  0(h6). 
The  formula  for  integration  is 

9 

f  f  F(x,y)  =  ^  ^  W1*Di  with  W2k  +  ]  =  25/324, W2k  =  10/81 

and  if  we  apply  this  to  the  integrals  for  X-  and  7  and  factor 'out  all  scale 
factors  we  get 

7  =  5  *  (Dl  -  D5)  +  4  *  (08  +  D2  -  D6  '-  D4) 

Y  =  5  *  (D7  -  03)  +  4  *  (D8  +  D6  -  02  -  D4) 

2  2 

and  use  abs  (X)  -  abs  (Y)  for  the  associated  magnitude.  If  we  sweep 
a  three-by-three  window  across  digitized  scene  D7  can  be  taken  as  the  upper 


10.  A.  H.  Stroud,  ibid. 


ORIGINAL  IMAGE  OF  DISK  (b)  VECTOR  FIELD  GENERATED  BY  (c)  EDGE  POINTS  OF  2(a)  FOUND 
WITH  A  ROOF  EDGE.  APPLYING  MOMENT  OPERATOR  BY  IDENTIFYING  THOSE 

TO  2(a).  POINTS  FOR  WHICH  /de=2  . 


left  hand  corner  while  D3  is  the  lower  right  hand  corner.  In  this  case  the 
direction  of  a  possible  edge  is  equal  to 

0  =  Arctan  (-■-——  +  tt/2 
Y  -  X 

b.  Calculation  of  the  Rotation.  The  vector  field  of  a  roof  edge  will 
look  like  the  vector  field  of  Figure  A-5.  To  find  roof  boundary  points,  we 
therefore  have  to  find  points  for  which,  in  a  neighborhood  of  such  a  point, 

/cd0  =  2tt 

The  smallest  region,  in  the  discrete  case  over  which  we  can  take  an  integral, 
is  a  two-by-two  window;  thus  our  algorithm  sweeps  a  two-by-two  window  across 
a  scene  and  computes  the  integral  /  d0  for  each  of  these  windows.  If  it 
turns  out  that  this  integral  is  equal  to  2n,  those  four  points  which  make  up 
the  window  are  classified  as  boundary  points.  To  calculate  the  integral  of 
a  two-by-two  window  we  use  an  approximation 

/d0  l  a£0 
l  =  1 

computed  by  a  computer  program.11 


For  the  purpose  of  this  experiment  the  procedure  used  to  generate  a  file 
which  is  the  file  of  detected  second  order  edges  in  the  following: 

1.  From  the  original  file  (scene)  two  files  are  generated;  one  (ACI) 
contains  SQRT  [(X)2  +  (Y)2];  and  the  other  (ANG)  the  angle  of  (0,  0<_  0<_  255) 
a  possible  edge. 

2.  From  the  ANG  and  ACI  files  one  new  file  AAA  is  created.  AAA  is 
created  by  sweeping  a  two-by-two  window  across  the  ANG  file.  The  rotation 
is  calculated,  and  if  a  point  is  classified  as  boundary,  then  to  the  corres¬ 
ponding  point  of  AAA  (initialized  at  zero)  Is  added  the  average  of  those  ele¬ 
ments  of  ACI  that  have  the  same  subscripts  as  those  of  the  two-by-two  window 
being  swept  across  ANG. 

Examples  of  how  this  method  works  are  shown  in  Figure  A-6. 


11.  R.  Machuca  and  A.  Gilbert,  "Finding  Edges  in  Noisy  Scenes,  IEEE 
Transactions  on  RAMI,  unpublished. 
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Section  4.  Evaluation 


The  methods  described  above  were  tested  on  disks  whose  edges  were  step, 
ramp  and  roof  edges.  The  step  and  ramp  edges  had  edge  height  equal  to  16 
while  the  roof  edge  was  constructed  by  beginning  at  the  center  with  gray 
value  equal  to  100  incrementing  by  one  to  gray  value  equal  132  and  then 
decrementing  by  one  to  gray  value  100.  All  files  were  128  x  128  x  8. 

To  test  the  effectiveness  of  the  different  operations  considered  here  we 
added  Gaussian  noise  of  different  standard  deviation  to  achieve  a  given 
signal  to  noise  ratio  and  then  tested  the  algorithms  (Figure  7). 

16  2 

The  SNR  ratio  was  measured  in  db;  that  is,  we  used  SNR  =  10  1°9-|q(^) 

n 

where  a  =  standard  deviation  of  the  noise.  For  the  ramp  and  step  edges 
we  usedmSNR  =  4,  5,  6,  .  .  .  ,  14  while  for  the  roof  edge  the  signal  to 
noise  ratios  used  were  10,  11,  12,  .  .  .  ,  20.  To  measure  the  effect¬ 
iveness  of  the  different  algorithms  we  graphed  PF  =  the  probability  of 
false  alarms  vs.  PD  =  the  probability  of  detection  (Figure  A-9).1 
Figure  A-8  contains  examples  of  processed  roof  edge  disks  with  SNR  =  13. 
The  graphs  of  PF  vs.  PD  (ROC  curves)  for  the  corresponding  operators 
appears  in  Figure  A-8. 

The  results  for  different  operators  and  step,  ramp  and  roof  edges  appear 
respectively  in  Figure  A-10.  These  graphs  show  that  the  performance  of 
the  moment  operator  is,  in  all  cases,  better  than  that  of  the  Sobel  opera 
tor.  A  significant  improvement  is  obtained  by  first  applying  the  average 
and  then  the  moment  operator.  When  the  signal-to-noise  ratio  is  high  the 
median  gives  better  results  than  the  average;  but  there  is  a  cross-over 
point  at  which  the  average  filter  gives  better  results  than  the  median. 


12.  X.  Abdou,  "Quantitative  Methods  of  Edge  Detection,"  Image  Pro¬ 
cessing  Institute,  University  of  Southern  California,  Los  Angeles,  Cali¬ 
fornia,  1978. 
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F i ^ 1 1 ro  A-p.  Poof  ednes  with  corresponding  ROC  curves. 
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